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Scattering on the energy shell is viewed here as the relation between 
the bound states of the Hamiltonian, restricted to sections on leads that 
are asymptotically independent, far away from the interaction region. The 
decomposition is achieved by sectioning this region and adding new leads, 
thus generating two new scatterers. So a resonant scatterer, whose S-matrix 
has sharp energy peaks, can be resolved into a pair of scatterers with smooth 
energy dependence. The resonant behaviour is concentrated in a spectral 
determinant obtained from a dissipative section map. The semiclassical limit 
of this theory coincides with the orbit resummation previously proposed by 
Georgeot and Prange. A numerical example for a semiseparable scatterer is 
investigated, revealing the accurate portrayal of the Wigner time delay by the 
spectral determinant. 

PACS numbers: 05.45. +b, 03.65.Sq, 03.80.+r 



I. INTRODUCTION 

Surfaces of section were conceived by Poincare as a means to decompose bounded classical 
motion into discrete maps of the section onto itself. Even if the longterm motion is chaotic in 
some sense, there are many important cases where the single return map is well behaved, so 
that the complex intermingling of trajectories results from repeated iterations. If a section 
is badly chosen, the first return map may well be singular and discontinuous with fractal 
boundaries between continuous regions. Evidently, such a Poincare map will not be a useful 
tool for the study of the classical motion. There are many systems for which no single 
surface exists that leads to a well behaved Poincare map, though these have not been much 
studied for obvious reasons. 

It is only with the work of Bogomolny [|TJ that surfaces of section where brought into 
quantum mechanics. At first this was viewed as an essentially semiclassical extrapolation 
of Poincare's method, while limiting the section to position space instead of the classical 
phase space. However, the approach of Rouvinez and Smilansky and Prosen [[|, allows 
us to define the transformation, undergone by the finite Hilbert space corresponding to the 
section, as the product of two unitary scattering matrices. These describe the two alternative 
unbounded motions, obtained by substituting either side of the system by a semi-infinite 
tube. Defining the scattering operators as Sl(E) and Sr(E), where E is the energy of the 
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open system, the condition for a bound state to exist at a certain energy is just that SlSr, 
or SrSl has a unit eigenvalue. 

The observation that a "bad" choice of section results in a maze of contributing orbits 
has a parallel in quantum mechanics that the correponding scattering exhibits multiple 
resonances. For example, consider the system in Fig. |l|(a). For the good choice of section, 
the classical orbits return to after a short time on the right side, whereas the return 
map for will be very complex. This corresponds to a smooth energy dependence of 
all the elements of the scattering matrix Sla (for the open system obtained by joining an 
open tube to the left of S^), as opposed to a spiky energy dependence of Slb due to many 
resonances. Let us now consider the case where the system under study is actually the open 
system depicted in Fig. |l|(b), obtained by joining a narrow open tube to the left of in 
system (a). The "natural section" that describes this new scattering problem would then 
be but having seen that this is a "bad choice" for the closed system, we would also like 
to bypass it in the description of the new open system (b). It may then be advantageous to 
decompose this open resonant system using a "good" section E^, just as in the case of the 
closed system. Indeed, we may consider the latter as the limit of a family of open systems 
in which the "bad" section is pushed evermore to the left until the opening at has zero 
width. The system (b) will then be decomposed into the auxiliary systems of Fig. [I|(c) and 
(d). The novelty with respect to the closed problem is that now there are modes entering 
the scattering region between and E# from the tubes both to the left and to the right 
of Fig. |l|(c). The important point is that the classical orbits remain only a short time 
in the scattering region, leading to simple mappings among the surfaces T*a and and 
nonresonant energy dependence of the S'-matrices. 

The purpose of this paper is to present the general method for the decomposition of a 
resonant scatterer into two less resonant scatterers by means of a surface of section. If the 
optimal decomposition renders both the resulting scatterers nonresonant, the objective is 
reached and the resonances will result from multiple scattering across the surface, as will be 
explicitly shown in Sec. |V]. If there is no choice of section for which one of the components 
is not resonant, we can clearly decompose it in the same manner and so on until all the 
components are nonresonant. A similar sequence of decompositions may be required for a 
bound system. There may not be any surface that avoids resonances in Sl or Sr. The 
further decomposition of either scatterer proceeds exactly as for initially open systems. 

First we review multichannel scattering theory on the energy shell, establishing the 
notation. In Sec. |TJ, we decompose the scatterer through a surface of section. The way 
that the resonant structure emerges from the decomposed scatterers is described in the 
following section. Particular attention is given to the time delay, as an overall indicator of 
resonance. Semiclassical approximations are discussed in Sec. [V[ These are not valid across 
a discontinuity, so that the local treatment for its traversal follows in Sec. |V|. We conclude 
with a numerical calculation for a semiseparable scatterer which allows us to evaluate the 
discontinuous approximation advanced in the previous section. We also verify that the 
localization of resonance peaks of the time delay is accurately determined by the determinant 
of a quantum Poincare map for the open system. 
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II. MULTICHANNEL SCATTERING 



We here review the theory of scattering on the energy shell so as to define the notation. 
Careful definition of the context determines the scope of the following theory in its present 
form. A general scatterer connects regions where asymptotic states, modes or channels, 
propagate independently. Without loss of generality, we may model the scatterer as a cavity 
connecting entrance leads ||. In many important cases, such as in mesoscopic devices, this 
structure is immediately evident, but even the scattering of waves from a localized object 
may be considered as radial backscattering along a tube of angular width 2tt. Defining xj 
as the longitudinal coordinate along the J-th lead and letting yj represent the transverse 
coordinates, the asymptotic region is specified by the requirement that, as xj — > oo, the 
Hamiltonian H(xj,yj,p x j,p y j) becomes independent of xj\ 

H(x j} yj,p x j,pyj) x ^° -^ L + Hj(yj,p y j) . (1) 

Here we assume that the scattering region lies near the origin of xj and that xj is always 
positive in the asymptotic region. We also assume bounded motion in the transverse direc- 
tions, implying that Hj(yj,p y j) has a discrete infinite spectrum of states \rij} with energies 

The Hilbert space for the stationary solutions of the Schrodinger equation are decom- 
posed into orthogonal subspaces in the asymptotic region of each lead. The J-th subspace 
is a superposition of states 

1 1 

{xjiVjM* = - 7 =exp(±ik nJ xj){yj\n J } = - 7 ={y J \n(x J )} ± , (2) 



where 



l2m(E - e nJ 



kn.j = \j h2 - ■ (3) 

In other words, the interaction between the transverse states \nj} is switched off as xj enters 
the asymptotic region. There remains only a trivial exponential dependence of the full state 
\rij) on the longitudinal coordinate. (We partially adopt the notation of Prosen in ||.) 

The states for which k n j is real are referred to as open channels, otherwise they are 
closed. For any given energy E, there will be a finite number, Aj, of open channels. The 
closed channels with rij > Aj take the form of evanescent modes, 

(xj,yj\nj) e = Ji exp(-\k nJ \xj){yj\nj} , (4) 



\J i\k n j 



for the wave to be bounded as xj — ► oo. However, we shall usually keep to (0) for both open 
and closed channels. 

A general solution \if>) of the stationary Schrodinger equation is decomposed in the J'th 
tube in terms of transverse amplitudes {nj\ipj} 
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(xj,yj\t/)) = —n= [{yj\ n { x j)} ++ { n Ai ) j} + {yj\n{xj)} {nj\tpj}] . (5) 

n=0 V fc nJ 

Here we could also have decomposed with the coefficients ± {n(xj)\ijj} to maintain the 
symmetry of the bracket notation, but it is easier to use a basis that is independent of xj. 
The incoming wave \ip{xj)}~ and the outgoing wave \ip(xj)} + at xj are thus defined as 

oo 

mx J )} ± = Y.\ n M} ±± i n Mj}, (6) 

n=0 

for the open channels. The generalized S matrix is now defined by the equation 

(^(xok, i^ 2 )} + , . . .) = s 0(*i)r , mx 2 )r, ...) . (7) 

Conservation of current among all the leads determines that the restriction of S to the finite 
block that includes all the open channels is a unitary matrix; time reversal symmetry implies 
that the full S matrix must be symmetric, i.e., Sy—Sji 0. Evidently, S depends on the 
choice of xj, but only up to a translation within the asymptotic region, implemented by a 
diagonal matrix in the \n(xj)} representation. The nontrivial energy dependence of S will 
be our main concern. 



III. DECOMPOSITION 



Let us consider an arbitrary scatterer, with J (> 1) leads such as sketched in Fig. 0(b). 
Cutting the scatterer by an arbitrary plane S , so as to avoid the leads, divides these into R- 
leads to the right of S and L-leads to the left (either of these sets may be empty). Likewise, 
we divide all the channels (open and closed) into i?-channels and L-channels. Henceforth it 
will be immaterial how the independent i?-channels (L-channels) are subdivided among the 
various i?-leads (L-leads). 

Define the coordinate Xq growing in the direction normal to S and positive on the i?-side 
and define the Hamiltonian Hq = H(xq = 0). Substituting H by Hq for all points to the 
right of So, we generate a semi-infinite lead. This can be joined onto H, to the left of So 
so as to form the L0 scatterer sketched in Fig 0(a). This procedure is reflected across S to 
form the OR scatterer, as shown in Fig. §(c). The L0 scatterer has only one 0-lead on the 
right and the states \n } are the eigenstates of H — p 2 xo /2m. These coincide with the states 
to the left of the OR scatterer. 

We now derive the relation of the L0 and the OR scattering to the original problem, i.e., 
the LR scattering decomposed by S . The constant profile of the H allows us to choose 
Xq — » in both the auxiliary scattering problems, which shall be assumed henceforth. The 
first step is to divide the scattering matrix for each of the scatterers into reflection and 
transmision blocks, generated by the subdivision of the the channels of each scatterer by S : 

a _ f r ll T lr \ _ f R° LL T L0 \ _ f R° RR T m \ , . 

\ J- RL -tlRR J \ 1 OL -n-00 / V 1 OR, -"-00 / 

To avoid confusion, we define the states related by each of these matrices 
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(\M*l)} + , = s (\Mxl)}-, \Mxr)} 

|^lo(^o)} + ) = S L (\Mxl)}~, \^lo(x )Y 
\*Pr(xr)} + , \*Por(xo)} + ) = S R (\ip R (x R )}~, \tpo R (xo)} 



(9) 
(10) 

(11) 

where R, L, OR, and L0 should be taken as specific instances of the index J in the definition 
(|) and in preceeding formulae. The matrices S, Sl and Sr account for the full scattering 
of three different systems. Our immediate task is to derive the elements of S from those of 
Sl and Sr. 

We have already identified the states on the right of the scatterer in Fig. 0(a) with those 
on the left of Fig. 0(c), as well as those on the left of Fig. 0(a) and (b). Now we match 
smoothly the wavefunctions of Fig. Q(a) and (c) at Xq = 0. Defining the operators 
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for any arbitrary power a, we have the matching conditions 

k 1/2 \ipLo(x )} + + k 1/2 \ip L0 (x )y = k 1/2 \iPor(x )} + + k 1/2 \il) 0R (x )Y 

and 

ikl /2 \ip L o(x )} + - ikl /2 \i/) L0 (x Q )}' = -ikjf 2 \iJ) 0R (x Q )} + + ikl /2 \il) 0R (x )y 
which determine uniquely 

|^oii(^o)} + = \^lo(x q )}~ 

and 

Inserting these equalities into equations (|IUJ) and flTTD, we obtain 

T 0L \^ L (x L )}- + R^ L0 (x )y = \^or(x )}- , 
T r\iPr{xr)}~ + R^\ip 0R (x )}~ = I^lo(^o)} - • 

This system of equations for l^/W^o)} - and \4>lo(xo)}~ is easily solved: 



fcOo)} - = [l - #£XoJ [RooTol\^l(x l )}- + T qr \i/; r {xr)Y 
IV^o)}" = [1 - -Roo^o]" 1 [RooTorIMxr)}' + T 0L \Mxl)Y 
where II is the unit matrix. Combining this result with (ID) and ([□]) again, leads to 

R° ll \Mxl)}- + 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 
(18) 



(19) 
(20) 



T L0 


1 - 


T3R r>L ' 
- XVqqJXqq 


-1 


'KTol\iI>l( 


x L )} 


+ T or \^r 




= \tpd 




ft-RR 


*pR 


>«)}- + 


-1 














Tro 


t - 


/ > L rfR 

- -Tt 00 -Tt 00 


RoqTor\*1>r 




' + T 0L \i/j l 




= \tpR 


(**)} 



(21) 
(22) 
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But this is just the required linear relation between the incoming waves 
\iPr(xr)}~ and the outgoing waves \4>l(xl)} + , \4'r( x r)} + , so that finally we have 



Rll 
Rrr 
Tlr 



R°LL + T 



R 



RR 



LO 



HO 



Jl 
1 



r>R r>L 

-^00 -"-00 

rtL rjR 

' -"do -"do 



-1 



r>R rp 



00 ± 0L , 



-1 



Rnr\T, 



00 1 OR 7 



T L0 
Trq 



r>R tdL 
- rl '00- rt 00 

rtL tdR 

- rt oo- rt oo 



-1 



OR 



0L 



(23) 
(24) 
(25) 
(26) 



In this way, the elements of the S matrix for the general scatterer are completely determined 
by those of the scatterers decomposed by S . 

We now discuss some elementary examples to illustrate the foregoing theory. The most 
trivial example of a scatterer is an arbitrary section of an open tube. Fixing the surfaces 
S# and as in Fig. |3|(a), we can "decompose" it with a surface So between these. The 
matrix Rqq = 0, whereas the transmission coefficients are diagonal: 



{ToR) nn 



exp(ik n d L0 ) 
exp(ik n d 0R ) 



K J- 0L) nn , 

(Tro)* ■ 



(27) 
(28) 



A major simplification results in the class of decompositions where there are no tubes 
on one side of S , as shown in Fig. ||(b). Then S reduces to 



LL 



R°LL + T 



L0 



a 



- rt oo- rt oo 



r*R rp 



'00 ± 0L 



(29) 



This will be the matrix studied in this Sec. [VI I|, since it already exhibits all the resonant 
structure of the general case. In the trivial case of a semi- infinite tube, Fig. 0(c), again 
R 



00 



and T 0L is given by (p7|), whereas 



(Roo) nn = -exp(2ik n d R ) 



(30) 



assuming Dirichlet boundary conditions at the end of the tube. 

It is important to remember that the generalized S-matrices defined in this section have 
infinite dimension. Only the blocks made up of all the open channels will be unitary. For 
this reason, the identification of the open channels for reflection with a unitary S-matrix 
in (^) is only valid because there are no open transmission channels. All projections of a 
unitary matrix into a subspace define a dissipative matrix with eigenvalues having moduli 
smaller than one, unless there is no interaction with the subspace that has been projected 
away 0. The addition of closed modes does not alter the dissipative nature of the reflection 
matrices in the general case where transmission is present. 



IV. RESONANCE STRUCTURE: MULTIPLE SCATTERING 

Let us assume that neither matrix Sr or Sl exhibits a sensitive energy dependence. As 
previously discussed, such a situation can always be reached in a sequence of decompositions. 
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(Even if the scatterer has a fractal structure, the wavelength settles the smallest scale of 
the fractal that can be resolved, and thus the number of auxiliary sections will be finite.) 
The sensitive energy dependence characteristic of a resonant scatterer then emerges from 
the possibility that one of the eigenvalues of Rq Rq or Rq Rq approaches the value one. 
In both cases the spectral determinant 

det(l - <X ) = det(l - i4< ) - . (31) 

As we show in || , any block with a smaller dimension than that of a full unitary matrix 
will generally have eigenvalues with moduli smaller than one. If the channels on the right 
are closed, Rq will be unitary, for the channels with real k®, but the product with R^ will 
still be dissipative. So we never have exactly zero in ( |3~TD for scattering situations; indeed 
these singularities characterize bound states. We may consider the bound states of closed 
systems as the limit of scattering systems in which the channels both to the right and to the 
left are closed. The condition (|3~1~D is then precisely the well established condition for the 
existence of an eigenstate 0-0]. Clearly, the order of the matrices Rq and Rq makes no 
difference to the bound state theory, though it does affect the full scattering problem. For 
the rest of this section this effect will be unimportant, so we shall abreviate both products 
Rq Rq and Rq Rq as simply R . 

The fact that Ro is a dissipative matrix allows us to expand 

oo 

(U - flo)- 1 = £ J*S . (32) 

We can thus interpret the formulae for each of the various scattering submatrices as the 
result of multiple reflections within the resonant region, each one contributing a term to the 
total scattering matrix. Note that ( |3"2D only converges for a scattering system. For a bound 
system, this is only a formal equality. Therefore, we can regularize the Bogomolny theory 
for a closed system by considering it as the limit of a scattering system. 

The fact that Rq is a well behaved operator with finite trace leads to the Fredholm 
expansion 

1 oo m 

(i-^^nEEW, (33) 

U m=0 fc=0 

where the spectral determinant 

oo 

D = det (11 - Ro) = ]T D m (34) 

m=0 

and 

i m 

D m = £An- fc Tri^. (35) 

This is a similar expansion to that used by Georgeot and Prange [[7], but their theory is 
a reworking of the semiclassical theory, as compared to the exact result. In practice, one 
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must work predominantly with the open channels contemplated by the semiclassical theory 
(see next section), but this finite matrix can be supplemented by as many closed channels 
as necessary to obtain convergence. 

The main point of the present theory is that we can define a resonant envelope for the 
energy dependence of all the matrix elements of S as the graph of D^ 1 (E). Certainly, there 
is a complex coupling among the different transverse states that varies smoothly with energy, 
so that this graph is not strictly an envelope. Even so, it is only at the peaks of D~ 1 (E) 
that we will find sharp resonant amplitudes of any element of the S matrix. Not all the 
peaks of D~ 1 (E) will manifest themselves for a single matrix element, but it is not surprising 
that collective properties of the S-matrix can be much more sensitive to this function of the 
surface of section. Indeed, the Wigner delay time, 

r(£) = -^^logdetS open , (36) 

where A is the number of open channels in S, measures the globally resonant nature of the 
scatterer, i.e., the peaks in r(E) correspond to the energies of longest permanence in the 
scattering region for some channel. Near such a peak, we may factor the rapidly varying 
part of S in (p9|), or in the more general CclSGS, clS 

S L ^D-\E)a L , (37) 

where er L has no peaks. It follows that near the resonant peaks 

r{E) » -^i og(jD -A det<T ) „ z^logD . (38) 

Thus, we obtain the peaks of t(E) « Tq(E), where we define the section time delay 

t (E) = -iftAi gdet(l - Ro)- 1 . (39) 

We shall verify the validity of this approximate form of the time delay for a numerical 
example in Sec. |V11| . 



V. SEMICLASSICAL APPROXIMATIONS 

The first step towards a semiclassical approximation of the preceeding theory is the 
truncation of the generalized S-matrices, which are restricted to the open channels. The 
reflections are then described by square sub-matrices, whereas the blocks accounting for 



transmission will be generally rectangular. Even though all the elements of (|23h26|) are 
then well defined, this approximation introduces spurious discontinuities in their energy 
dependence as each new channel is opened. The semiclassical theory can only avoid these 
by some form of complex continuation of the classical motion. It should be noted that, 
even so, the discontinuities in energy do reflect some of the qualitative features of the full 
quantum theory, so that energy integrals may be quite reasonable. In any case, we expect 
that truncation will be a good approximation in an energy range with a constant number of 
open channels, in the limit when this number is large. 
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We now associate the transverse eigenstates \rij} singled out by EBK-quantization (see 
e.g., H) on a given section Sj to the invariant tori for the classical Hamiltonian Hj(yj,p yj ) 
defined by ([I]). In other words, we consider the canonical transformation {yj,p yj ) — > (Ij, 9j), 
such that Hj = Hj(Ij), obtained from the multivalued generating function F v (Ij,yj): 



dF u 



dF u 



Pyj 



e. 



dyj dlj 

(where the index v distinguishes the different branches). Then we approximate 

d 2 F„ 



(40) 



{yj\nj} « 



dljdyj 



exp [iF v (Ij,yj)/h] 



with 



Ij{n) = ft(n + /x/4) 



(41) 



(42) 



where /i is the Maslov index for the torus. (The phase differences on passing to a different 
branch are additive, having been included in F v .) The number of open channels Aj(E) is 
the largest integer n satisfying the condition 



Hj{Ij{n)) < E . 



(43) 



Let us first study scattering where all the points in Ej, subject to (£3), return to the 
same region. Semiclassically we have a reflection, S = R, corresponding to a conservative 
classical map. Each quantized torus determines a tube of trajectories that re-intersect Sj, 
i.e., it "evolves" while preserving its shown in Fig. |j. The overlap of the evolved 

torus state with the basis of the quantized tori is just 



{n'j\R\nj} « 



d 2 F v 



dljdl'j 



expiiF^ijj'jyn] 



(44) 



where the canonical transformation is generated by 

dF u nl dF v 



dl'j 



7 J 5 



dlj 



0, 



(45) 



This general formalism for semiclassical maps was originally developed by Miller [H. The 
branches, indexed by u, correspond to each transverse intersection of the new torus with the 
basis torus Ij = h(n + ///4). Both the new torus and the basis tori are closed curves (for 
nonresonant scattering) so there will be an even number of branches. Caustics, leading to 
spurious singularities in (fHf), occur when the intersection of the tori is non-transverse, i.e., 
where two intersections coalesce. In this region the matrix element should be expressed in 
terms of Airy functions instead of (§4|) (see e.g., |§). As we shall discuss below, this simple 
picture of a returning torus that overlaps with the basis of tori is fragmented into a fractal 
mosaic for a very resonant system. However, we can now keep to a simple description for 
each nonresonant scatterer into which the section decomposes the original resonant system. 
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It is important to note that we need not worry about the problem of quantizing a compact 
phase space. Our torus basis is infinite, though we are concerned with the projection within 
a classically invariant region (which grows with E). The matrix elements among the open 
channels in the semiclassical approximation depend only on the classical dynamics within 
this region. The discrete action variables form a very privileged basis. We cannot transform 
R to the (semiclassical) \8j} basis without doing a Fourier sum over elements {n'j\R\nj} 
that lie outside the open region, not defined semiclassically. The same difficulty involves 
passing to the \t/j} representation of R. However, the semiclassical evaluation of unitary 
transformations relies entirely on the stationary phase approximation. We can thus make 
ordinary semiclassical changes of bases within the allowed region by assuming that there are 
no stationary phase points outside it. 

For scattering where the open channels are not restricted to a single lead, we must 
consider the conservative classical map defined on the union of several surfaces of section, as 
shown in Fig. [|. Again, the classical motion takes place on the region ([|3J) for each section, 
so that the allowed region will vary from section to section. The basis of eigenstates for each 
section corresponds to a set of quantized invariant tori of that particular section. 

To determine the semiclassical approximation to the ^-matrix when there is more than 
one lead, we again consider the tube of trajectories defined by one of the quantized tori in 
one of the sections. But now this tube will generally split, reintersecting the various sections 
along open segments. Semiclassically, there will be nonzero matrix elements with all the 
quantized tori with which the segments intersect. The elements of the reflection block are 
still given by ([44]), whereas the transmission elements from the J- lead to the i^-lead depend 
of the generating functions F u (Ij,Ik) for the canonical transformation (Ij,6j) — > (Ij,6j) 
in the form 

1/2 

exp[iF„(Ij,I' K )/h] . (46) 

By adopting fixed action-angle variables along Sj, with the normal coordinate xj, the 
action for orbits returning to the same lead, 

F u (Ij, I'j) = £ 9jdlj + f Pxj dxj , (47) 

except for Maslov indices, is evaluated along the classical orbit joining the torus Ij to I'j. 
This can also be obtained by taking 

Fu{yj, y'j) = J Pyjdyj + j p XJ dxj (48) 

and then evaluating the change of basis by stationary phase. To evaluate the transmission 
actions, we first note that the difference between fl44|) and (|46"D is only that Ij and I' K belong 
to different torus bases, rather than they belong to different sections. Defining F° as the 
generating function for the canonical transformation corresponding to this change of basis 
we obtain 

F v {I Jt I' K ) = F°{Ij, I' K ) + [ K 9 K dI K + [ Px dx , (49) 

JIk(Ij) J'Ej 



{n' K \T\nj} 



d 2 F u 



dIjdI' K 
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where 



OF 

at = ~ e « < 50 > 

describes the original torus with action variable Ij in the (Ik, 6k) coordinates; Ik(Ij) is the 
action variable of the orbit on this torus that will arrive at I' K and the variable x in the last 
integral is assumed normal to both sections. 

So far, this review of semiclassical scattering theory applies to an arbitrary scatterer. The 
problem arises, for a very resonant scatterer that the classical map fragments into arbitrarily 
small regions which cannot be quantized when their area is smaller that Planck's constant 
@. However, we have seen that it is always possible to decompose a resonant scatterer by 
appropriate surfaces of section. It is then possible to adopt semiclassical approximations for 



each of the components and then to evaluate the blocks of the full S matrix given by (|23| - p6D 
in the stationary phase approximation. 

For a typical resonant scatterer, cut by a "good" section, E , the classically allowed 
area in So will be considerably larger than the union of the i?-leads and of the L-leads 
projected onto So- The classical map corresponding to the various T blocks of the component 
scatterers will be a simple injection of orbits into So, or its time reverse. If we decouple this 
part of the S matrix, there results the view of semiclassical scattering as essentially that 
of the So surface onto itself with projections onto entrance and exit regions, as previously 
proposed ||. The semiclassical neglect of any backscattering in the LO or RO scattering 
depends on the smoothness of the Hamiltonian. It is certainly inadequate for a discontinuity 
in the Hamiltonian within the range of the longitudinal wavelength. This case will be treated 
in Sec. |VI[ 



If we apply the Fredholm theory in Sec. |T| to [1 - RqqRqq]' 1 and to [1 - R^R^' 1 



nioJ 



and then use the semiclassical expressions for each of the reflection matrices, we succeed in 
rederiving the scattering theory of Georgeot and Prange ]7| from first principles. This is 
especially valuable, because it was originally obtained by the rearrangement of the semiclas- 
sical contributions to the scattering matrix, which is not valid in the resonant context. The 
spectral determinant (133) will be described by the classical periodic orbits that cross the 



surface S rearranged into composite orbits or pseudo-orbits. As we have seen, the spectral 
determinant is never zero, because the R matrices are dissipative. No matter how long lived 
the typical classical orbits may be, as counted by the number of traversals of S , there will 
be a cuttoff for the period of the periodic orbits included in the Fredholm theory. This is 
limited to the dimension of the R matrices. 

To conclude this section, we note that in the important case of systems with hard walls, 
i.e., open billiards connected by leads as in Fig. |^, the quantized tori for any section, central 
to the foregoing theory, will be defined as the level curves of Py. The phase change between 
the two branches of will depend upon the choice of Dirichlet or Neumann boundary 

conditions. It may well be advantageous to use the perimeter of the closed billiard (£# in 
Fig. []), corresponding to the classical Birkhoff map (or bounce map) instead of a section 
such as So- This was the point of view adopted in ||. In this case, all scattering will be 
considered as a reflection across However, in this case, it is essential to take into account 
the discontinuity on entering the billiard, which is the subject of the next section. 
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VI. DISCONTINUITIES 



If the Hamiltonian is discontinuous across a given plane S, we should construct two 
sections E^ and E# immediately on either side of E, as shown in Fig. This is an important 
example of a scattering system which is clarified by the use of extra surfaces of section. In 
the following theory we only address the local problem of scattering from E^ to T,r. With 
respect to Fig. |7], we may consider the I — > R scattering as decomposed by E L , or that E R 
decomposes the L — > r scattering. Thus the full I — > r scatterer will be decomposed in two 
stages. In the simple semiseparable example discussed in Sec. [VI \ both the I <-> L and 
R «-> r propagators are trivial. 

To solve the scattering problem across the discontinuity between the sections E l and E# 
involves the same smoothness conditions as (|T3J) and , except that now the wave vectors 
fci? 7^ k l and we have different basis states {tir} and \ul} to match on either side. Recalling 
the definiton of the operators k" in (fHJ), we have 



k 

ik 



-1/2 



1/2 



k 

ik 



-1/2 



1/2 
R 



\M + + HrY 



(51) 
(52) 



Decomposing the S matrix into reflexion and transmission blocks (H), then yields 



k 



-1/2 



k 



1/2 



(H + Rll)\iI>l}~ + Tlr\iJ>rY 



-1/2 
R 

1/2 
R 



[SL + RrrMr}- 

-tl + RrrMr}- 



Tr L \^ l Y 
•Trl^lY 



(53) 
(54) 



The fact that these equations are valid for any \4>r}~ or \ipi}~ allows us to equate 
separately the operators acting on either of these functions. However, to derive matrix 
equations, we must transform between the natural bases on either side. Thus, defining the 
orthogonal matrix 



(55) 



we obtain 



*Z ,/! (i 
fe! /2 (n 



Rll) 
Rll) 



Uk-R 1/2 T RL 
UkfTRL ; 



(56) 
(57) 



where the matrix Rll is in the L-representation, whereas (T^) ran / 
natively, we obtain from ([54]) and the transpose of (|55| ) that 



{n R \T RL \n' L }. Alter- 



k R 1/2 (l + R 



k 



1/2 
R 



RR) 

(1 - Rrr) 



LR 



U T kl l/2 T 
U T k L /2 T LR 



(58) 
(59) 



where, again, the indices L,R specify the bases. Notice that the set of equations (^,^) (and 
its solution) is related to the set (|56|j57| ) by interchanging L,U <->• R, L/" T (time reversal). So 
we restrict ourselves to (|58| , p9D , which can be solved for 
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Tlr = 2 (ktfu* k~ L 1/2 + k R 1/2 U T k 1 / 2 )' 1 , (60) 
R RR = § (kTu T k- L ^ - k- R ^U T kT) T LR ; (61) 

with T R l and Rll given by time-reversal as we previously noted: 

T RL = 2{kTuk R ^ + kl^Ukfy\ (62) 
Rll = \ (k L /2 Uk R 1/2 - k- L 1/2 Ukf) T RL . (63) 

It is straightforward to verify that the S matrix for the discontinuity is symmetric. The 
verification of unitarity is also simple, but lengthy ||. 

The matrix inversions in the preceding formulae will be calculated by truncating the 
bases on the right and on the left of the discontinuity A possible choice is to limit U 
to the open channels on either side, but their numbers may not be the same. A better 
alternative is to notice that the matrix elements U n r n decay exponentially when there is no 
intersection between the quantized tori corresponding to \til} and to \n' R }. Therefore, the 
open channels should be supplemented with states which intersect with them on either side 
of the discontinuity. 

It is interesting that the preceeding argument relies on a semiclassical criterion of torus 
overlap to calculate U n i n = {n' L \n R }, even though we use this to calculate classically for- 
bidden reflections. We shall use a similar semiclassical approach to estimate the channels 



exhibiting resonances in Sec. |VII| . This semiclassical inspiration will now be pushed fur- 
ther to suggest an approximation that avoids the matrix inversions in S. To this end, we 
notice that the semiclassical evaluation of the nondecaying U n / n depends on the points of 
intersection for the corresponding tori (such as shown in Fig. [T0| ). We could then view 
each intersection as the projection of an orbit transmitted or reflected at the discontinuity. 
Thus each orbit contributing to a given matrix element faces the same discontinuity in the 
longitudinal Hamiltonian, as depicted in Fig. |8]. What is the probability amplitude for a 
wave \n R } incident on the right to propagate to the left channel \n' L }l Separating out the 
longitudinal motion, we have 

\x,y R ) = k~ R /2 exp(-ik nR x) + (r RR ) n , n k~ R /2 exp(ik nR x) {y\n R } , x > , 

\x,y L ) = (t LR ) n > n k^£ 2 exp(-ik n > L x) {y\n' L } ,x<Q. (64) 

Here, tiR and r RR are the "longitudinal" transmission and reflection matrices. By requiring 
that both waves match at the discontinuity at x = we obtain: 



2 1 </£Nfe' ' (65) 





W*= hF-ir 1 • (66) 



which are the well known results for the one- dimensional step potential. In the case of a 
two-dimensional discontinuity we still have to take into account the coupling U n > n between 



13 



transverse modes at each side of the discontinuity. The expression for the approximate 
transmission matrix T is thus 



LRJn'n 



LR)n'nU n 'n 



(67) 



Similarly {r RR ) n ' n U n ' n gives the probability amplitude for not being transmitted to chan- 
nel n' . As this probability is distributed over the L basis we must switch back to the R 
representation. Our approximate result for the reflection matrix reads 



Rrr = U T [r RR U] . 



(68) 



The corresponding approximation for the transmission and reflection matrices for waves 
inciding from the left are obtained by the time reversal operation: 



Rll 



t RL U T 



U 



vllU 1 



(69) 
(70) 



We can clarify the nature of this two dimensional version of a "sudden approximation" Q8"7D 
by noting that 



k~ L 1/2 fk R /2 U T 



k]l 2 Tk^^U J 



-1/2 TT T 



(71) 



as compared to (|6"0"D, so that T amounts to a reordering of the operators that define T. This 
allows the simplification of not having to invert the matrices in the exact formula flBTf). 



VII. A SEMISEPARABLE NUMERICAL EXAMPLE 

We now study a simple though nontrivial example of the decomposition of a resonant 
scatterer. This is chosen as the semiseparable system sketched in Fig|| i.e., the potential 
V(x,y) is separable for x > and for x < 0, but discontinuous at x = 0: 

f \uj\y 2 x < 
V(x,y) = l \uly 2 U<x<a . (72) 
[ oo x = a 

We choose ojl > u R , so that the classical motion is broader on the right than on the left: 
the outline in Fig. |9] represents an equipotential curve. (The main difference with respect to 
the system studied by Prosen 0] is that it is open in the left.) Computations were carried 
out with the values u L =l, 0^=0.71423, a=30, h—1, and mass m—1 (in appropriate units). 

Separability on the left allows us to bring the entrance section to x = 0~. The 
position of the Poincare section is also arbitrary, because of separability, so we bring So to 
x = + . There being a single open lead, the decomposition of the scattering matrix is given 
by (p9|), where the reflection matrix on the right of So is diagonal, with elements 

{n R \Ii£ \n R } = - exp(2ik Rn a) . (73) 
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The transmission matrices Tlo and T l are just those for the passage through a dis- 
continuity of the potential discussed in the last section. This permits us to evaluate the 
accuracy of the approximation proposed there within the resonant theory. The reflection 
Rq may also be calculated in this approximation. In spite of its motivation in terms of 
classical orbits, one should notice that this "sudden approximation" mixes the motion that 
fails to escape before it is propagated diagonally by Rq . Thus, the full Poincare reflection 
matrix Rq Rq is not diagonal, in contrast to the classical motion which is integrable until 
it escapes. 

The S matrix was calculated in an energy range corresponding to 15 open channels, but 
excluding an interval of 0.2 near the thresholds of the 15th and 16th channels. Evidently 
each channel threshold is given by 

E n = hu L {n - §) . (74) 

The number of transverse states on the right corresponding to open channels grows from 21 
to 22 in the energy window, with our choice of the frequency ujr. 

We can predict which entrance channels may exhibit resonance behaviour by reverting 
to the qualitative semiclassical picture of Sec. [V], even though the discontinuity prevents an 
accurate semiclassical calculation. The modes in the open lead correspond to the concentric 
circles of in Fig. |H] with area lisin^ — \), whereas the states in the cavity correspond to the 
ellipses with area 2tx{tir — |) and semiaxes with the ratio ojr. Semiclassically, an incoming 
mode only excites the inner states corresponding to ellipses that intersects its circle. If 
these ellipses, in their turn, only intersect circles with < 15, there should be straight 
nonresonant backscattering for this mode. It is easy to see that, by this criterion, it is only 
for til > 15cj|, « 7 that the resonant structures should arise. 

In Fig. 11 we compose D = det[l — Rq Rq ] and the section time delay with Wigner's 
time delay. It is impressive how the overall structure of this global imprint of the resonances 
is re-composed from both nonresonant matrices Rq and Rq . It is important to note that 
the peaks in the time delay are only poorly correlated to the bound states of the separable 
system obtained by closing the lead at x — 0, also shown in Fig. [II]. [A comment about this 
plot is in order. Recalling that in (^) we had discarded the slowly varying part of the time 
delay, we shifted t q (E) by a fixed amount so that both average times (r ) and (r) coincide. 
Then logarithms were taken and, finally, the curves corresponding to t q (E) and D(E) were 
shifted downwards to make comparisons easier.] 



In Fig.[L2| we show the resonance structure of individual S matrix elements. Evidently, 
each one of these does not exhibit all the resonance peaks in the time delay, but the de- 
terminant D does limit the positions allowed for these peaks in all cases. We find that the 
positions of the resonances are reasonably obtained, except for very fine structures, but the 
amplitudes of the individual matrix elements are not predicted by the spectral determinant. 
In these figures we also compare the sudden approximation with the exact calculation. This 
approximation works remarkably well for both nonresonant and resonant channels up to 
i=13. 
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VIII. CONCLUDING REMARKS 



The spiky energy dependence of the S matrix characteristic of resonant scattering in 
quantum mechanics corresponds to complex orbital structure of the classical limit. In both 
theories these complications can be explained in terms of the multiple iterations of a rela- 
tively simple mapping defined on an appropriate surface of section. The quantum theory 
for such an open system is based on the pioneering developments of Bogomolny, Prosen 
and several papers by Smilansky and coworkers on bound systems, but here we have the 
advantage that the map is dissipative. Thus, our formula (|33D converges and we could, in 
principle, regularize the section theory for bound systems by considering them as the limit 
of a family of open systems. 

By focusing the scattering problem on a section map, there emerges the central role of 
the spectral determinant. As the openings of the scatterer are closed, the complex zeroes 
of D(E) converge onto the real eigenergies of the resulting bound system. We have shown 
that even for the open system the full profile of this energy function is found to portray the 
qualitative features of the scattering time-delay. Our numerical example reveals that, even 
for individual elements of the S matrix, the possible resonance peaks are restricted to those 
of the spectral determinant. 

The picture of the on-the-shell scattering as relating sections on various leads, which may 
be decomposed by splicing the scatterer and adding more leads, has a clear semiclassical 
interpretation. For a scattering system with two degrees of freedom, there is only one 
freedon left transverse to the leads. The corresponding classical Hamiltonian is therefore 
integrable and this allows us to associate an invariant torus (closed curve) of the asymptotic 
Hamiltonian to each scattering channel. The classical propagation of each torus through the 
scatterer generally breaks up these tori and the elements of the S matrix result from the 
intersection of their fragments with the exit channels. This theory goes back to Miller, but 
its resummation based on a section was achieved by Georgeot and Prange. The advantage 
our new derivation of this theory is that it proceeds from first principles, rather than as a 
rearrangement of Miller's orbit sum. Only in this way can we show that the approximate 
semiclassical spectral determinant is based on a dissipative rather than a unitary matrix. 
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FIG. 1. A single return of the orbit to corresponds to multiple iterations of the Poincare map 
for E^. The surface of section decomposes the resonant scatterer in (b) into two nonresonant 
scattering systems (c) and (d). Accordingly, the scattering trajectories are complex in (b) but 
simple in both (c) and (d). Thick curves represent a single contour of the potential energy. 



(a) 



L 




FIG. 2. The multilead scattering system. The section So decomposes the full scatterer (b) into 
the nonresonant systems (a) and (c). Thick curves represent a single contour of the potential 
energy 
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(a) 



'LO 



'OR 



*L Z *R 




(c) 



'LO 



'R 



FIG. 3. Three scattering systems: a straight open tube (a), a cavity connected to a single lead 
(b), and a tube closed by a hard wall (c). 





FIG. 4. Single-lead scattering viewed on the surface of section. The circles represent channel 
tori selected by EBK quantization. If scattering is nonresonant a torus (dotted) is mapped into a 
closed curve [thick line in (a)]. In a resonant case (b), the mapping is discontinuous and the torus 
returns to the section as a set of fragments. This set may be infinite and even fractal. 





FIG. 5. In multi-lead scattering the full surface of section is the union of the sections of all leads 
(two in this example). Shown is a torus in the left section [dotted, (a)] which is divided into a 
reflected part and a transmitted part [thick lines in (a) and (b), respectively]. 
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FIG. 6. A billiard connected to leads. Plot are a "standard" section So and the Birkhoff section 
Tjb- Thick lines represent the contour of the open billiard. 




FIG. 7. A scattering system discontinuous across a surface E. The sections E^ and Y,r allow 
for isolating the problem of the discontinuity in the scattering from E/ to E r . 









K n'L 




k 2 
K nR 








> 



E 



I 



FIG. 8. The "equivalent" one dimensional scattering problem for the discontinuity. A particle 
with energy E propagates from the right channel nR to the left channel n'L through the disconti- 
nuity E. The jump in the longitudinal momentum is k n >L — k n R, where Ti 2 k^, L /2m = E — e n >L and 
h 2 k 2 nR /2m = E- e nR . 
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a 



-o 



FIG. 9. The semiseparable system consists of two waveguides: a semi-infinite one (x < 0) and 
a finite one (x > 0). Dirichlet conditions are assumed at x = a. Also shown are the auxiliary 
sections and Sq. 




n R =10 



> y 



FIG. 10. Quantized tori on each side of the semiseparable system. Circles correspond to channel 
modes and ellipses to cavity modes. The outermost circle corresponds to the last open channel. 
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FIG. 11. Logio-linear plot of the exact time delay t(E), the section time delay tq(E), and the 
spectral determinant D(E). The spikes at the bottom represent the energy levels of the bound 
system obtained by closing the cavity with a hard wall at x = 0. 
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14.6 14.8 15 15.2 15.4 14.6 14.8 15 15.2 15.4 



E 

FIG. 12. Elastic cross sections \Su\ 2 as a function of energy. Full lines correspond to the exact 
calculations. Dots represent the results of the "sudden approximation" (see text). 
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